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ABSTRACT 

We investigate the far-infrared properties of X-ray sources detected in the Chandra Deep 
Field-South (CDF-S) survey using the uhra-deep 70 (a.m and 24 [im Spitzer observations taken 
in this field. Since only 30 (i.e., « 10%) of the 266 X-ray sources in the region of the 70 \xm 
observations are detected at 70 \xm, we rely on stacking analyses of the 70 (a.m data to char- 
acterise the average 70 i_im properties of the X-ray sources as a function of redshift. X-ray 
luminosity and X-ray absorption. Using Spitzer-IRS data of the 5w//r-BAT sample of z « 
active galactic nuclei (hereafter, AGNs), we show that the 70/24 ^m flux ratio can distin- 
guish between AGN-dominated and starburst-dominated systems out to z w 1.5. Among the 
X-ray sources detected at 70 \xm we note a large scatter in the observed 70/24 [im flux ra- 
tios, spanning almost a factor of 10 at similar redshifts, irrespective of object classification, 
suggesting a range of AGN:starburst ratios. From stacking analyses we find that the average 
observed 70/24 [im flux ratios of AGNs out to an average redshift of 1.5 are similar to z « 
AGNs with similar X-ray luminosities {Lx = lO'*^^'*^ erg s^^) and absorbing column densities 
(Nh ^ 10^^ cm^^). Furthermore, both high redshift and z « AGNs follow the same tendency 
toward warmer 70/24 |^m colours with increasing X-ray luminosity (Lx). From analyses of 
the Swift-BAil sample of z ~ AGNs, we note that the 70 |xm flux can be used to determine 
the infrared (8-1000 |i.m) luminosities of high redshift AGNs. We use this information to show 
that Lx= lO'*^^'*^ erg s^' AGNs at high redshifts (z = 1 — 2) have infrared to X-ray luminosity 
ratios (hereafter, Lir/Lx) that are, on average, 4.7^2V ^^'^ ^■^■'^-2 6 t™ss higher than AGNs 
with similar X-ray luminosities at z = 0.5 — 1 and z ~ 0, respectively. By comparison, we 
find that the Lir/Lx ratios of Lx= 10'^^-'^* erg s"' AGNs remain largely unchanged across 
this same redshift interval. We explore the consequences that these results may have on the 
identification of distant, potentially Compton-thick AGNs using Ljr/Lx ratios. In addition, 
we discuss possible scenarios to account for the observed increase in the Ljr/Lx ratio with 
redshift, including changes in the dust covering factor of AGNs and/or the star formation rates 
of their host galaxies. Finally, we show how deep observations to be undertaken by the Her- 
schel Space Observatory will enable us to discriminate between these proposed scenarios and 
also identify Compton-thick AGNs at high redshifts. 

Key words: galaxies: Active Galaxies, High Redshift, Infrared: Galaxies, X-rays: Galaxies 



1 INTRODUCTION 

Deep X-ray surveys undertaken by the Chandra and XMM- 
Newton observatories have provided the most efficient method to 
date of identifying large numbers of active galactic nu clei (here- 
after, AGNs) out to high redshifts (z ~ 5; e.g.. Alexander et"al] 
I2OOII : iBarger et al.l l2003l : iBauer et al.l |2004| : ISzokolv etal1l2004 
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iBrandt & Hasingeil |2003 : iMainieri et afl |2005|) . Due to the high 
penetrating power of hard (i.e., >2 keV) X-rays, they provide a 
means of identifying AGNs that is much less biased toward un- 
obscured AGNs than, for example, that obtained at optical wave- 
lengths alone. Consequently, deep X-ray surveys have revealed that 
a large fraction (up to ~ 75%; Mainieri et al .l2002l : lDwellv & Page! 
[2OO6; Tozzi et al. 20oA lTaier et alj2007h of highlv obscured AGNs 
are often uniden tified or misclassified by observations at other 
wavelengths (e.g. iBarger et aI.ll2003l : ISzokolv et alj2004h . Further- 
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more, as the central engines of AGNs are the only compact ob- 
jects capable of producing luminous X-ray emission (i.e., Lx > 
10"*^ erg s^'), observations at these high energies provide the most 
reliable method of measuring the intrinsic power of the central 
engine without suffering from significant contamination by star- 
formation. 

Although X-rays provide an ideal means of studying the 
central engines of AGNs, they yield limited information on the 
regions outside the central few parsecs. For example, although 
X-ray observations have been used extensively to measure the 
amount of absorbing gas (i.e. A'h) along ou r line of sight to the 
central engin e (e.g. Mushotzky et al. 1 9931 ; ISassani et al] 1 19991 : 
iMalizia et al]|2003l : iGuainazzi et alJIlOOSn Tthev tell us little about 
the spatial extent and covering factor of the putative obscur- 
ing dusty torus, a key component of the unification models that 
have become a cornerstone of AGN research over the past two 
decades (e.g.|An tonucci & Miller 1985; Krolik & Begelman 1988; 
I AntonuccH 1 1 993h . Furthermore, as any X-ray emission from stel- 
lar processes will be swamped by the presence of an AGN, X- 
rays alone provide little insight into the relationship between star- 
formation and AGN activity that is predicted by both g alaxy for- 



mation models (e.g. .Kauffmann & Haeh nelt 20 00; Be n son et al 
20031; iGranato et all 120041 ; iBower et alj l2006riBooth & Schave 



20090 and locally defined black hole-bulge relationships (e. 



Magorrian etalll 19981; iGebhardt et al.ll2000l ; iTremaine et alj|2002l ; 
McLure & DunloDl l2002Y Observations at infrared (hereafter, IR) 



wavelengths can provide significant insights on both these counts. 

Like X-rays, the longer IR wavelengths (i.e., >5 \im) are 
largely unaffected by absorption from interstellar gas/dust. How- 
ever, because IR radiation is typically produced via the reprocess- 
ing of shorter wavelength light by dust, it also provides a means 
to study the dusty environments surrounding AGNs. Furthermore, 
as both AGNs and starbursting systems are capable of producing 
large amounts of IR radiation, the study of AGNs at these wave- 
lengths can provide insights into the processes that connects AGN 
activity and star-formation. However, despite considerable gains in 
our understanding of the IR properties of AGNs since the launch of 
the Spitzer Space Telescope, there remains many unresolved ques- 
tions. For example, we currently lack a clear picture of how the 
IR emission of AGNs is affected by the differences in the intrinsic 
power of the central engine and, perhaps more importantly, whether 
it evolves over cosmic time. As a result, our understanding of both 
the putative torus and the interplay between AGN activity and star- 
formation is limited. 

To date, most studies of the IR properties of AGNs have con- 
centrated on the relatively short wavelength data provided by the 
IRAC instrument (3.6 ^.m, 4.5 \Lm, 5.8 |J.m and 8.0 |xm) and the 
24 [im observations taken wit h the MIPS instrument on-board the 
Spitzer Space Tele s cope (e.g. |Lutz et al. | |2004 IStemetalJ|200^ ; 
iDonlev et al]|2007l ; lOaddi et al.ll2007t) . However, studies of high 
redshift objects at these shorter wavelengths are hindered by the 
presence of strong spectral features at A,A.< 12 ^.m (in particular 
poly cyclic aromatic hydrocarbons [PAHs], silicate absorption fea- 
tures and strong emission lines) that shift into these wavebands at 
z > 0.7. For example, the IR luminosities of high redshift galaxies 
are overestimated by a factor of ~ 2 — 5 (depending on redshift and 
Ljr) when using 24 yim data alone compared to that calculated us- 
ing the full IR SE P incorporating longer wavelength data (at 70 ^m 
and 160 ^.m; see lPapovich et al j|2007l ; iMagnelli et al.ll2009l Mag- 
nelli et al. 2009 [in prep.]). Furthermore, as emission from high red- 
shift starbursting systems peaks at longer wavelengths than those 
probed by the 24 ^m waveband, the sole use of near and mid- 



IR observations hampers our ability to efficiently distinguish be- 
tween star-formation and AGN dominated systems. This too can 
be resolved by incorporating longer wavelength IR data into our 
analysis. Undoubtedly, the launch of the Herschel Space Tele- 
scope, with its ability to observe the universe at 55 — 672 [im, will 
help in addressing these issues. However, progress can be made 
now by applying stacking analyses to deep Spitzer 24 [im and 
70 ii.mobservations, enabling us to measure the average IR prop- 
erties of AGNs at a sensitivity level comparable to that achievable 
for individual sources in theproposed deep Herschel fields (e.g. the 
HGOODS Key ProgrammfQ; PI: D. Elbaz). 

In this paper we investigate the mid to far-IR properties of 
X-ray detected AGNs that lie in the confusion-limited 70 (xm 
and 24 \im Spitzer-MIPS Far-Infrared Deep Extragalactic Legacy 
Survey (FIDEL) and Great Observatories Origins Deep Survey 
(GOODS). We identify AGNs using the 1 Ms Chandra Deep Field- 
South X-ray observatio ns (which have hi gh quality X-ray spectral 
constraints published in lTozzi et al ][2006h and explore their 70 \xm 
fluxes, luminosities, and 70/24 fim flux ratios (hereafter, 57o/524) 
as a function of redshift. X-ray luminosity, and X-ray absorption. 
As part of our analysis we use the Swift-BKT sample of local, well 
studied AGNs to develop diagnostics to distinguish between star- 
burst (hereafter, SB) and AGN dominated systems at z < 1.5. In 
the CDF-S field we find evidence that more X-ray luminous AGNs 
have warmer, more AGN dominated IR colours and that high-z, 
low X-ray luminosity (i.e., z = 1 — 2, Lx = lO^^^^"* erg s^') AGNs 
are significantly more IR luminous than their low redshift counter- 
parts (i.e., z = 0.5 — 1, Lx = lO^^^^"' erg s^'). This increase in IR 
luminosity may be interpreted as tentative evidence of larger dust 
covering factors in high redshift AGNs, supporting results from 
deep X-ray surveys which find evidence of larger lev els of A'h at 
similar redshifts (e.g. iLa Franca et al]|2005l ; lHasingei ii2008. How- 
ever, a proportion of this increase may also be attributed to higher 
levels of star-formation in these systems. We show that deep, far- 
IR Herschel surveys will resolve the dominant process behind this 
increase in IR luminosity and enable the discrimination between 
AGN and star-formation dominated systems out to z ~ 5-6. 



Throughout this work we adopt Hq = 11 km s Mpc 



-I 



= 0.27, and Slj^ = 0.73. We report lo errors for detected 
sources and provide 3o upper limits for undetected sources (<3o). 



2 DATA 

Our main focus throughout this work is the analysis of the mid 
and far-infrared (MIR [5-30 \xm\ and FIR [30-300 vim], respec- 
tively) emission of X-ray detected sources in the CDF-S field. To 
aid in the interpretation of the CDF-S data, we also make use of 
archival Spitzer-[RS spectra and IRAS data for the Swi/Z-BAT sam- 
ple of local X-ray AGNs, two archetypal AGNs (NGC 1068 and 
NGC 6240), and a distant quasi-stellar object (QSO) sample. We 
describe the CDF-S sample in SjZT]and the comparison samples in 



2.1 CDF-S Sample of X-ray Detected AGNs & Starbursts 

2.1.1 Chandra Data 

The X-ray data for our m ain sample are taken from the 1 Ms 
CDF-S X-ray observations jGiacconi et alj|2002h . as analysed by 

' URL: |http://herschel.esac.esa.int/Docs/BCPOT/GOODS-HerscheLp^ 
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Figure 1. Rest-frame hard X-ray luminosity (i-2-lOkeV-' rsdshift (z) of the AGNs and SBs identified in the 1 Ms CDF-S X-ray observations, separated 
according to column density (Nh). Those sources detected at 70 [am are indicated with large circles. As described in i|2.1.1l this sample contains biases in 
favour of more luminous sources and higher column densities at higher redshifts. To mitigate the effects of these biases we define a more heterogeneous 
subsample (the 'Restricted' sample), contained within the grey shaded region. Coloured histograms indicate the number of AGNs in each redshift and X-ray 
luminosity bin that have column densities Nu < 10^^ cm^^ (green) and A'h = ICf^^'^ cm"^ (red). The black histogram shows the distribution of the full 
AGN-l-SB sample. Solid colours/lines are used to indicate the restricted sample, hatching/dashes indicate the (unrestricted, flux limited) parent sample. The 



curved dotted line indicates the constant, unabsorbed hard X-ray flux (6.8 x 10 
X-ray AGNs. 



erg s cm ) that best fits the Lx-z distribution of the 70 |jm detected 



lAlexander etai] j2003h . In total there are 201, 304, and 326 X-ray 
detected objects in this region down to limiting 2-8 keV fluxes of 
10"'^, 10"^^ and IQ-^'' erg s"' cm^^, respectively. The CDF-S 
data is deep enough to detect star-formation activity out to z ~ 1 . 
Since the primary aim of this study is to exp lore the IR proper- 
ties of X-ray detected AGNs, we have used the lBauer et alj |2004|) 
X-ray source classifications to separate AGNs from star-forming 
galaxies. In the full catalogue of X-ray sources, there are 288 X- 
ray AGNs and 15 X-ray SBs (the remaining 23 are classed as either 
normal galaxies or stars and are excluded from further discussion). 
The absorption corrected 2-10 keV X-ray luminositi es (Lx) and 
X-ray absorption column densities (Nn) are taken from lTozzi et al.l 
(2006). Here, we only consider those 266 X-ray AGNs/SBs with 
well defined Lx and A'h measurements. Spectroscopic and photo- 



metri c redshifts for the se 266 X-ray AGNs/SBs were also taken 
from lTozzi et al.l J2006l) : 137 hav e spectroscopic re dshifts (113 of 
which are described as secure, see iTozzi et alj2"oo3 for details). We 
use photometric redshifts for the remaining 129 X-ray AGNs/SBs 
that lack spectroscopic redshifts. 

In Fig. [T] we present the Lx-redshift distribution of the 266 X- 
ray AGNs/SBs considered in this study. As is to be expected in flux 
limited samples, there is a strong bias toward the detection of more 
luminous sources at higher redshifts. Furthermore, as is shown in 
Fig. (2] there is a bias toward higher column densities at higher X- 
ray luminosities because a) only the brightest X-ray sources can be 
detected behind A'h > 10^^^^^ cm^^and b) it is difficult to measure 
low values of A'h at high redshifts due to the absorption cut-off be- 
ing shifted out of the lowest Chandra energy band. In our analysis 
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Figure 2. Rest frame hard X-ray luminosity (i2-10keV) versus column 
density (Nh) of AGNs and SBs identified in the 1 Ms CDF-S X-ray ob- 
servations. Also shown is the Lx-Nn distribution of the Swi^-BAT sample 
of z AGNs which cover the same region of the parameter space as the 
CDF-S sample. The selection criteria for inclusion in the 'Restricted' sam- 
ple is indicated by the shaded region. 



we consider two main samples: a) the full sample of classified X- 
ray AGNs/SBs and b) a 'Restricted' sample limited to those AGNs 
withLx = lO'*^"'^'* erg s"' andA'n < lO^"' cm"^ within the redshift 
range of 0.5 < z < 2.0 (indicated by the shaded regions in Figs[T] 
and[2l(. The 'Restricted' sample is used to mitigate the effects of 
absorption and selection biases on the average results; the Lx and 
Nil range of the restricted sample is also well matched to the z ~ 
Swift-BAT comparison sample (see i]2.2.1| >. 



2.1.2 Spitzer Data 

The 24 ^m and 70 ^m 5p;rzer-MIPS GOODS/CDF-:fl observa- 
tions (PID: 20147; RL: D. Frayer) and FIDElB legacy surveys 
(PID: 30948; PL: M. Dickinson) are the deepest available at these 
wavelengths and therefore present the best datasets to study the 
MIR to FIR emission of AGNs in the high-z universe. W hile the FI- 
DEL region covers the full E-CDFS jLehmer et alj2005l) . the deep- 
est region is the 10' x 10' region centred on the CDF-S field cov- 
ered by the GOODS observations. The raw 70 \im data were pro- 
cessed off-line using the Germanium Reprocessing Tools (GeRT), 
following the techniques described in lFraver et al.l J2006l) . The fi- 
nal mosaics, produced by MOPEX0, have a pixel scale of 4.0". The 
70 \im 3a sensitivity is 2.0 mjy in the CDF-S field and ~3.0 mjy 
in the outer E-CDFS regions. The 24 |j.m images produced using 
MOPEX have a pixel scale of 1.2". The sky coverage at 24 \im. is 
inhomogeneous and the sensitivity ranges from 30 to 70 (xJy (5ct). 
Assuming a SNR of 4.1 for the the 24 \im observations (corre- 



^ URL: http://www.stsci.edu/science/goods/ 
3 URL: http://ii-sa.ipac.caltech.edu/data/SPITZER/ FIDEL/| 



sponding to an mean flux limit of ; 
tected sources in the FIDEL field. 



^50 ^Jy) there are 14025 de- 



URL : |ht^j^ssc£pitzer£altech£d^ 



2.1.3 X-ray-IR source matching 

We matched the classified X-ray sources to the 70 \xm sources as- 
suming a search radius of 2.0" for those sources in the GOODS 
field with 24 [xm/IRAC counterparts and 4.0" in the larger FIDEL 
field. This yields 30 matches which are listed in Table [T] repre- 
senting « 11% of the 266 X-ra y AGNs/SBs. Using the P-statistic 
{P = \ — exp(— Jt«6^); e.g., see lDownes et al.ll 19861 and Eqn 1 in 
IPope et ani2006l) and a 70 \im source density n ^l.Ox lO-' deg^^, 
the chance of one or more 70 \im source lying within 4.0" of an X- 
ray source is P ~ 0.2%. Given the small value of P, it was not nec- 
es sary to apply the corr ection to the matching procedure described 
bv lDownes et alj h986h . 

We matched the X-ray sources to the 24 \im sources using a 
3.0" search radius; this search radius was found to provide the best 
compromise between finding real matches while reducing the num- 
ber of spurious matches. Out of the 266 CDF-S X-ray AGN/SBs we 
found 172 have 24 \mi counterparts (i.e., ~ 65%). The number of 
spurious matches in the CDF-S field was estimated to be ~ 4.5% 
using two approaches ( 1 ) by calculating the area of the field covered 
by 24 fi.m sources, assuming that each source is a circle of radius 
3.0", and (2) by shifting the positions of the 24 \x.m sources by a 
random displacement between 20" and 55" and then re-matching 
to the X-ray sample. This spurious matching fraction is consis- 
tent with the P-statistic (~3%), assuming n ~1.5x 10^ deg^^ (see 
above). In spite of the larger P-statistic associated with the 24 \xm 
m atches compared to t he 70 |i.m matches, the correction described 
in lDownes et al. I ( [198^ need only be applied to 15 (6%) of the 266 
X-ray sources and has negligible effect on any of our results. 

All but one of the 70 vim detected sources has a 24 \xm coun- 
terpart. Further investigation reveals that the lack of a matched 
24 ^im counterpart is the result of the larger PSF of the 70 |i.m 
image; at the corresponding position in the 24 \xm image there is 
a cluster of sources that lie just outside the 3.0" matching radius 
which blend to form a single "source" in the 70 ^m image. To mit- 
igate source blending we de-blend the 70 ^.m sources with Gaus- 
sians placed at the 24 ^m positions. It is these 'de-blended' fluxes 
that are reported in Table [T] 

Finally, one of the 70 \im detected X-r ay sources lies within 
3.0" of two 24 \im sources jTozzi et al.l2006i Index: 3 1); in Table[T] 
we list both 24 \xm matches but assume that the brighter of the two 
24 \xm sources is the real match in our analysis. 



2.1.4 24 \im and 70 \im stacking procedure 

As the majority of the X-ray sources are not detected at 70 \im 
we rely on stacking analyses to provide insight into their average 
MIR to FIR proper ties. Stacking was performed using the code of 
iHuvnh et al.l ( l2007b . Cutouts of 128" x 128" (32 x 32 pixels) were 
made from the 70 \im image at the stack positions and combined 
using a weighted mean. The 70 nm stacked flux density was mea- 
sured using an aperture of 8.0" at the stacked image center. Off- 
set stacks were generated using random but nearby offset positions 
(< 64") in the 70 nm image. Two hundred randomly offset stacks 
were made and the uncertainty in the stacked flux is taken as the 
standard deviation of the measured flux density in the 200 offset 
stacks. 

We created stacks that either include or exclude the 70 \im. 
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detected sources; the former (hereafter, "all") provides a global 
average of the FIR properties of all X-ray sources and the latter 
(hereafter, "undetected") provides the average of the 70 \im. unde- 
tected X-ray sources that is less skewed by individual bright 70 |i.m 
sources. We also stack only the 70 \im. detected sources for com- 
parison (hereafter, "detected"). Because the all stacks often appear 
to be dominated by the few 70 ^.m detected sources we generally 
use the detected and undetected stacks in our analysis. 

In order to isolate trends between the IR emission from 
AGNs and other physical properties we further split the 'Re- 
stricted' sample (see [jZTTJ in terms of Lx (= 10'*^^^"' erg s^', 
= 1043-44 gj-g ^-1^^ 10^2 cm-2, = 1022-23cm-2) and red- 

shift (z = 0.5 - 1.0, z = 1.0 - 2.0). 

2.2 Comparison Samples 

In our analysis we compare the results from the sample of CDF-S 
X-ray AGNs/SBs to the IR properties of two arche typal AGNs, a 
sample of bright QSOs from lRichards et sA\ |200^ and, most im- 
portantly, a sample derived from local (z < 0. 1) AGNs in the Swift- 
BA T catalogue of har d X-ray detected AGNs (see i]2.2.1| appendix 
and lTueller et"ai]|2008,) . For the latter we selected those AGNs that 
cover the same range of X-ray properties (i.e. Lx and A'h; see 
Fig. |2} as the CDF-S X-ray AGNs sample and also have archival 
Spz/zer-IRS spectra and/or IRAS flux density measurements. The 
result is a sample that is directly comparable to the more distant 
CDF-S AGNs, but which has far superior IR data due to the rela- 
tive proximity of the sources. 

2.2.1 The Swift-BAT Sample 

The Swift telescope is currently undertaking a survey of X-ray 
bright (Fi4_i95 jj^ev > 5x 10^" erg s^' cm^^) AGNs using its on- 
board Burst Alert Telescope (BAT); the secon d data release (DR2 ) 
of this survey has recently been published in iTueller et al.l j2008l) . 
As the BAT instrument is sensitive only to very hard X-rays (14- 
195 keV) this survey has resulted in a sample of locaj^ AGNs that 
is largely unaffected by absorption (to A'h < lO^'* cm^^) xjjg Lx 
and A'h dis tributions of the 104 AGNs for which a rchival 2-10 keV 
data exist ( iBassani et Jill 9991 ; Iwiiiter et al.ll2009t) are comparable 
to the CDF-S sample (see Figs. 7 and 10 of lWinteret alj|2009l) . 
Of these 104 AGNs, 61 have flux density measurements from all 
four IRAS bands (these 61 objects are hereafter referred to as the 
BAT//i?A5 sa mple) which we convert to Ljr using Eqns. 2 and 3 
in Table 1 of I Sanders & Mirabe] ( Il996l) . A search of the Spitzer 
archives reveals that 36 of the 104 AGNs with 2-10 keV data have 
5p;rzer-IRS spectroscopy between 5.2-38 vim (these 36 objects are 
hereafter referred to as the BAT/IRS sample) which we use in con- 
junction with IRAS data to derive the expected observed 24 \im and 
70 [i-m fluxes at the redshifts covered by the CDF-S sample. All of 
the AGNs in the BAT/IRS sample are also in the BAT/IRAS sam- 
ple. See appendix for full details of the \RSIIRAS analysis. 

On separating the BAT/IRS sample into SB and AGN domi- 
nated objects we find that the Sjq/S2'\ ratio is able to discriminate 
between these two types of sources out to z ~ 1.5 (see shaded re- 
gions in Fig.|4]and appendix). Furthermore, we note that the 70 |i.m 
flux density can predict the infrared luminosities of all the BAT/IRS 
objects (i.e., both SB-dominated and AGN-dominated AGNs), as 

^ 135 of the 153 AGNs wit h measured redshifts in the BAT sample (pre- 
sented in lTueller et alj2008l) have z < 0. 1 . 



well as the starburst galaxies described in lBrandl et al. I ( l2006l) , to 
within a factor of ~3. This compares to a factor ~12 uncertainty in 
LjR using 24 \im flux densities alone (see appendix and Fig. |A3b . 

In the plots that follow, we indicate the average IR properties 
of the BAT/IRS sample using red and blue lines for the AGN and 
SB dominated systems, respectively. Shading is used to indicate the 
range of MIR and FIR properties of the BAT/IRS sample. 

2.2.2 A'GC 1068. NGC 6240 and the QSO sample 

In addition to the BAT/IRS sample we also compare the CDF-S 
AGNs with two well studied heavily obscured AGNs, NGC 1068 
and NGC 6240. The former is regarded as the "quintessential" type- 
II AGN, showing evidence of a hidden broad line region in po- 
larised light (e.g. Antonucci & Miller 1985); the latter is a heavily 
obscured AGN that is SB-dominated at infrared wavelengths and is 
often cited to characterise the properties of faint X-ray AGNs. The 
tracks for these AGNs were derived from archival Spitzer-\KS spec- 
tra following the same procedure as used for the BAT/IRS sample 

as outlined in the appendix. 

Finally, we use the average QSO SED of iRichards et al.l 

( l2006l) to give an indication of the typical IR properties of lumi- 
nous, unobscured type-I AGNs. 

We stress that the comparison samples considered here do not cover 
the full range of SEDs seen in the general galaxy population. It is 
known, for example, that very strongly starbursting systems have 
larger Siq/Sia flux ratios than NGC 6240 (e.g. Arp 220). Fur- 
thermore, some quiescent galaxies are known to have low 570 /S2A 
flux ratios similar to those of AGNs in the BAT/IRS sample (see 
Alexander et al., in prep). However, our main aim here is to ex- 
plore the IR emission of AGNs with X-ray properties comparable 
to the CDF-S AGNs, which excludes such extreme systems. We 
find that known starbursting systems, such as those presented in 
lBrandletal.l ( [2006) . all produce expected 570 /S2A flux ratios either 
consistent with or higher than those of the SB-dominated BAT/IRS 
AGNs, as indicated by the shaded region in Fig. lA3l 



3 RESULTS 

In our analysis we investigate the IR properties of the CDF-S X- 
ray AGNs and SBs, including the S-jq/Sia flux ratios, the IR lu- 
minosities and the IR— X-ray luminosity ratios. To aid in our anal- 
ysis we compare the properties of the CDF-S X-ray sources with 
expectations based on the more detailed analysis of nearby AGNs 
with the same X-ray properties as the CDF-S sample (i.e., Lx; A'h). 
Since only a minority of CDF-S X-ray AGNs are detected at 70 \im, 
our analyses rely significantly on 70 \im stacking analyses and the 
24 \xm fluxes of individual sources. 

3.1 Individual 70 (xm detected sources 

3.1.1 X-Ray properties and redshift distribution 

In Fig. [T] we show the X-ray luminosity-redshift distribution for 
the CDF-S X-ray AGNs/SBs and highlight the 70 ^m detected ob- 
jects. As might be expected, there is a preference for bright X-ray 
AGNs to be detected at 70 [im. However, there are up to 2 or- 
ders of magnitude difference in Lx for 70 tim-detected sources at 
the same redshift, suggesting a large range of X-ray/IR luminos- 
ity ratios. As we show in Fig. [S] detection at 70 [xm is strongly 
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comparison AGNs, normalised to Lir=10'^ Lq, at z = 0.25 — 2.5. 



dependent on the 24 \im flux density. None of the AGNs with 
S24 < 200 (i-Jy, and fewer than a quarter (10/45; ~22%) of those 
with 200 < 524 /^Jy< 500, are detected at 70 vim. Conversely, the 
majority (19/28; Ri68%) of those AGNs with 524 > 500 ^Jy are 
detected at 70 i^m. 



The 70 \im detected sources span a broad range of X-ray 
properties. Of the 29 sources detected at both 70 \im. and 24 ii.m, 
15 (Ri52%) are non-broad-line AGNs (non-BLAGNs), 6 {^2\%) 
are broad-line AGNs (BLAGNs) and 8 (fa28%) are X-ray detected 
SBs. These 70 vim detected sources comprise 15/218 (~7%) of the 
non-BLAGNs, 6/34 (Ril8%) of BLAGNs and 6/14 (Ri43%) of the 
X-ray detected SBs; the bias toward the BLAGNs and X-ray de- 
tected SBs could be due to the BLAGNs and SBs having high Lx 
and large IR-X-ray luminosity ratios, respectively. The 70 |i.m de- 
tection of AGNs shows no apparent dependency on absorbing col- 
umn density: for example, at z < 1.5, 6/17 35%) 70 vim detected 
sources are unabsorbed (A'h < 10^^ cm^^), 5/17 (~ 29%) are ab- 
sorbed (Wh ~ 10^2 - 10^3 cm"2) 5/^7 35%) heavily 
absorbed (A'h > lO^-' cm"^). 



3. 1.2 Distinguishing between SB and AGN dominated systems 

In Fig. |4] we show the S-h)/S2A flux ratios of all of the CDF-S 
X-ray AGNs/SBs and compare them to the range expected for 
the BAT/IRS AGN sample and the well studied local AGNs and 
quasars; see section 12.2.11 and appendix. The S-10/S2A flux ratio 
provides an effective discrimination between those AGNs with 
AGN-dominated or SB-dominated IR SEDs to z 1 - 1.5. Us- 
ing the tracks derived for the BAT/IRS sample we find that 15/23 
(^65%) and 5/23 (Ri22%) of the 70 ^m detected CDF-S AGNs 
have 570/524 flux ratios consistent with SB -dominated and AGN- 
dominated systems, respectively; the remaining three 70 vim de- 
tected AGNs lie in regions of the 570/524— z plot that are consis- 
tent with either SB or AGN dominated systems. Ten (~ 42%) of 
the 70 vim detected CDF-S AGNs have Sjo/Sil flux ratios within 
a factor of 1.5 of the NGC 6240 track, five of which are X-ray 
absorbed and five of which are X-ray unabsorbed. The maximum 
570/524 flux ratio of the 70 vim detected sources is 24. 4 ± 1.9. This 
is marginally inconsistent with iPapovich et alj ( l2007h who found 
3/30 (fa 10%) X-ray detected AGNs in their 70 vim selected sam- 
ple with 570/524 > 30. However, the larger area covered by their 
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survey increases the likelihood of finding rare AGNs with more 
extreme 570/524 flux ratios that are not represented by our compar- 
ison samples; see ii2.2l 

Among the 70 \im detected sources we find no difference in 
S70/524 flux ratio between BLAGNs and non-BLAGNs, nor be- 
tween X-ray absorbed and X-ray unabsorbed AGNs (i.e., A'h ^ 
10^^ cm^^ and A'h < 10^^ cm^^, respectively), implying that the 
material that absorbs photons at X-ray and optical wavelengths is 
optically thin to rest frame MIR-FIR radiation, see i]3.1.1| The av- 
erage 570/524 flux ratios and redshifts of the absorbed (unabsorbed) 



70 \xm detected AGNs are 57o/524 = 8.2±4.3 (11.8±7.3) and 
z= 1.1 ±0.6 (0.76 ±0.560 respectively. 

All six 70 \im detected X-ray SBs have S10/S2A flux ratios 
consistent with that expected of SB-dominated systems. However, 
one of the SBs has a 57o/524 flux ratio larger than that expected 
from the BAT/IRS sample and is typical of the extreme SB systems 
presented in Brandl et al. (2006), see Fig. lABI The 70 \im detected 
X-ray SBs also span the same range of Sj[)/S24 flux ratios as the 
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seven 70 yim detected X-ray AGNs that lie in the same redshift 
interval (i.e., z = 0.2 — 0.7), providing further evidence that these 
AGNs have SB-dominated IR SEDs. 



3.1.3 IR fluxes and luminosities of individual, 70 \im-detected 
AGNs 

Since we can provide a good characterisation of the basic IR SEDs 
of the 70 nm-detected CDF-S X-ray AGNs, we can accurately es- 
timate their IR luminosities. We derive Ljr for the 70 nm detected 
sources by first selecting one of the five tracks plotted in Fig. |4] 
that best matches the Sjo/S24 flux ratio of a given CDF-S X-ray 
AGN, then scaling its IR luminosity to reproduce 570 of the CDF-S 
source if observed at the same redshift. Estimating Lir in this man- 
ner assumes that the IR SEDs of AGNs at the redshifts covered by 
the CDF-S observations are similar to those seen in the comparison 
AGN samples. We take some confidence that this is indeed the case 
as the range of 570 / ^24 flux ratios of the 70 [im detected subsample 
is largely bounded by the tracks of our comparison sample. 

In the left panel of Fig. [5] we show the Lir of the 70 ^.m de- 
tected sample plotted as a function of Lx- We find no significant 
difference in Ljr between BLAGNs and non-BLAGNs or between 
X-ray absorbed and unabsorbed AGNs. This independence of Ljr 
on classification provides further supporting evidence that the ma- 
terial that absorbs photons at X-ray and optical wavelengths is op- 
tically thin to rest frame MIR-FIR radiation, see ^ j3TT] & [3TIl 
The 70 \Lm detected CDF-S X-ray AGNs are typically more IR 
luminous than their counterparts in the BAT/IRAS sample of local 



AGN. This result may be due to the low sensitivity of the 70 [im ob- 
servations, which could lead to a bias towards the detection of the 
most IR-luminous AGNs; indeed, the range of IR luminosities (~ 2 
orders of magnitude) for the CDF-S X-ray AGNs is narrower than 
the range of X-ray luminosities (~ 3 orders of magnitude). Over- 
all, the Lir/Lx ratios of the 70 \im detected CDF-S AGNs/SBs 
span almost 2 orders of magnitude. 

The 70 nm detected X-ray SBs have higher Ljr/Lx ratios 
than the majority (20/22; ^ 91%) of the 70 ^.m detected X-ray 
AGNs. We find all three 70 ^.m detected X-ray AGNs with Lx < 
10'*' erg s^' have Lir/Lx ratios and Sjg /Sy^ flux ra tios compa- 
rable to the 70 \xm detected X-ray SBs jTozzi et alj l2006 index; 
525, 538, 575; see Table[Tll. However, two of these s ources (index; 
525, 5 38) have optical spectra consistent with AGNs jSzokolv et al.l 
|2004 while the remainder (index; 575) has a flat X-ray spectral in - 
dex in the deeper 2 Ms CDF-S catalogue (F < 0. 1; ILuo et alj2008l) . 
While this confirms their classification as AGNs, it is clear that their 
IR SEDs are SB dominated. 

As only a small fraction of the CDF-S X-ray AGN sample 
is detected at 70 [im we are unable to reliably constrain the IR 
properties of the majority of the CDF-S X-ray AGNs/SBs using 
the 70 \im data alone. To better characterise the IR properties for 
the majority of the CDF-S X-ray AGNs we therefore use stacking 
analyses. 

3.2 24 \im and 70 stacking analysis 

The advantage of stacking analyses is that it provides the average 
properties of sources that lie below the individual source detection 
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limit. However, naturally, information is lost on the fluxes of in- 
dividual sources. We also stacked the 70 \im and 24 (i.m data fol- 
lowing the procedure outlined in i]2. 1.41 and present the results in 
Tabled although the majority of the CDF-S X -ray sources are de- 
tected at 24 ^j.m, we stacked the 24 [im data to provide flux density 
constraints consistent with those obtained at 70 vim. For complete- 
ness, we report the results from stacking the whole sample as well 
as a range of subsamples (split into bins of Lx, A'h- 24 |j.m flux 
density, and object classification). As noted in i]2.1.4| we stacked 
all X-ray sources in each subsample ("all"), and also stacked only 
the 70 \im detected objects ("detected") and only the 70 [im un- 
detected objects ("undetected"). The majority of our analyses will 
focus on the 'Restricted' sample to reduce selection and sensitivity 
biases (see fjZTTTJ. 

All of the CDF-S X-ray sample is detected significantly in the 
70 |i.m stacks, with an average flux density of 57o=4824 ± 210 (xJy. 
Significant detections (S/N > 3) are also found from 70 [im stack- 
ing analyses for the majority of the subsamples; however, we note 
that only 9 (~ 50%) of the eighteen 70 urn undetected subsamples 
are significantly detected. From stacking the CDF-S X-ray sample 
in bins of 24 \i.m. flux density (i.e. ^24=50-150 jxJy, 150-500 (xJy 
and ^500 viJy), we find that the 70 |i.m flux density and the signif- 
icance of the detection is positively correlated with the 24 \xm flux 
density of the stacked sources, confirming our previous result for 
the 70 \im detected sources; see i]3.1l 

We find that when the 'Restricted' sample is stacked in bins of 
Nn and z (i.e., stacks 5, 6, 7 and 8 in Table|2} only one 70 |J.m unde- 
tected stack is detected at 70 \un at a significance > 3cT. The upper 
limits for these non-detections reveal no significant differences in 
the IR properties of these subsamples; they are all consistent with 
S7o/524~4 and LirR^ lO" Lq. The lack of any significant differ- 
ence between these stacks may be interpreted as providing further 
tentative evidence that the material that absorbs X-rays is optically 
thin to rest frame MIR~FIR radiation, see ^ ^3TT][3X2l & [3T3l 
Similarly, none of the Lx = 10"*^^^"* erg s^', 70 \im undetected 
stacks (rows 13 and 14 of Table |2l( are detected at > 3cT in either 
redshift bin (z = 0.5-1 and z = 1-2). However, in this case the up- 
per limits do provide significant insight into the IR properties of the 
CDF-S X-ray AGNs, revealing that high Lx AGNs have, on aver- 
age, warmer S-1Q/S24 flux ratios than their low Lx counterparts (as 
discussed in i]3.2.1b . 

We proceed with the analysis of the stacked data following 
a similar procedure as outlined in i]3.1l for the 70 \im detected 
sources. 

3.2.1 SB/AGN contribution to the average IR emission of AGNs 

In the lower panel of Fig.|4]we show the 570/524 flux ratios of the 
'Restricted' subsamples split in terms of Lx and z (rows 11, 12, 13 
and 14 of Table|2](. We find that the average S70/524 flux ratios of 
all four 70 |i.m undetected stacks lie within the range of S70/S24 
flux ratios expected for the BAT/IRS AGNs out to z ~ 2.5, and run 
roughly parallel to the average AGN-dominated and average SB- 
dominated tracks (see i]2.2.1| and appendix). The simplest interpre- 
tation of this result is that there is little change in the average IR 
colour with redshift for AGNs in both the Lx = 10'*^^^"' erg s^' 
and Lx = jq'*-'^'*^ erg s^' bins. However, the conservative up- 
per limits on the S-jq/Sia flux ratios of the Lx = lO'*-'^'*'* erg s^' 
70 ^.m undetected AGNs place them below those of the Lx = 
JQ42-43 gj-g ^jj^ undetected AGNs in each redshift bin. 

Based on the BAT/IRS tracks this result suggests that the more X- 
ray luminous AGNs may have, on average, more AGN-dominated 



IR SEDs than their lower Lx counterparts; however, deeper data 
will be required to confirm this result. However, we find that the 
stacks of the Lx = 10^^^^-^ erg s^' AGNs lie in ambiguous re- 
gions of the 570/524— z plot, which limits the conclusions that can 
be directly derived on the relative contributions from AGN and SB 
activity to these stacks using the S-JQ/S24 flux ratios alone; in Sj4T| 
we explore a variety of approaches to constrain the relative AGN 
and SB contributions using additional data. 

3.2.2 Average IR fluxes and luminosities 

The average IR luminosities of the stacked subsamples are calcu- 
lated by taking the same approach as that used for the individual 
70 \im detected sources (see i]3.1.3b . For those stacks with less than 
3CJ detections we use the nominal 570 /524 flux ratio (rather than 
upper limits) when determining the closest match out of the five 
tracks considered; although we note that using upper limits to se- 
lect the appropriate track changes the estimates of Lir by less than 
20%. In both redshift bins the nominal 570/524 flux ratios of the 
Lx= lO'*^"'*^ erg s"' and Lx= lO"^^""*"^ erg s"' 70 ^m undetected 
stacks are most closely matched by the average AGN and average 
QSO SEDs, respectively. We note, however, that the 3 <J upper lim- 
its on the 570/524 flux ratios of stacked Lx= lO"*^^^"* erg s^' AGNs 
are also consistent with the average AGN SED; see Fig.|4] 

In Fig. |5] we plot the IR luminosity versus Lx for the CDF- 
S X-ray sources and the BAT//iM5 sample. In the case of the 
Lx = lO^^^^"' erg s^' undetected AGNs, we find that Ljr is larger 
by a factor of ~ 5 for the z= 1.0 — 2.0 AGNs when compared to 
the z = 0.5 - 1.0 AGNs (Lir =2.5 xlO" Lq and 5.0xI0'° Lg, 
respectively), and a factor of ~20 times larger than the average 
42 <log(Lx/erg s^'X 43, %< lO^^ cm^^ BAT/IRS AGN (i.e., at 
z ~ 0; Lir =1.3x lO"^ Lg,). This is in spite of these subsamples hav- 
ing almost the same average X-ray luminosities (Lx =4.9x10^^, 
4.3 X 10"^^ erg s^'and 3.7 x 10^2 erg s^'for z = 1.0-2.0, z = 0.5 - 
1.0 and z ~ AGNs, respectively). This difference corresponds to a 
factor of 4.T^^Q^ and 12.7^2 6 increase in Lir/Lx from z = 0.5 — 1 
to z = 1—2 and from z~Otoz = 1.0 — 2.0, respectively (these con- 
servative errors correspond to the full range of IR luminosities de- 
rived using all five tracks in Fig.|4j. The stacked data therefore sug- 
gest that the average Lir/Lx ratio of Lx = lO'*^^'*-' erg s^' AGNs 
strongly evolves with redshift, as illustrated in Fig. [6] Although this 
result is based on stacking analyses, we show in i]3.3| that we obtain 
the same result using 24 ^.m constraints for individual sources. 

The 70 vim undetected stacks of the Lx = jq'*-'^'*^ erg s^' 
AGNs at both z = 0.5-1 and z = 1-2 are not significantly detected, 
limiting the constraints that we can place on Lir using the 70 \x.m 
data alone; the limits are Lir <1.5x10" Lq and <4.6xl0" Lq 
for the AGNs at z = 0.5-1 and z = 1-2, respectively (these upper 
limits include the uncertainty in the 570-LiR correction). However, 
because the upper limits on the 570 /524 flux ratios of these stacks 
provide significant constraints to the range of potential matching 
BAT/IRS AGN SEDs we can estimate the average IR luminosi- 
ties from the stacked 24 ^im flux densities to a much higher de- 
gree of accuracy than would normally be the case (to within a fac- 
tor of ~3, rather than ^^12). This reduced uncertainty is based on 
the fact that the 570/524 flux ratios of the Lx = 10'*^"'*'* erg s"' 
AGN stacks restricts us to AGN-dominated BAT/IRS SEDs, which 
all have similar 524-Lir ratios. Using the stacked 24 \im flux 
density we obtain Lir=5.6x lO'*' Lq and 1.2xl0" Lq for the 
Lx = 10^3-44 ergs-1 AGNs in the z = 0.5 - 1.0 and z = 1.0- 2.0 
redshift bins, respectively, corresponding to a factor of ~2 in- 
crease between these redshift ranges. For comparison, the aver- 
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with iLutz et al.l ( |2004|) (i.e., ~ 3 and ~ 2, respectively^ To cal- 
culate the vLv(6/vm)/Lx ratios of the CDF-S AGNs we assume 
the same SEDs as those used to derive their Ljr/Lx ratios (i.e., 
Average AGN and Average QSO for Lx = 10^^^"*^ erg s^' and 
Lx = lO^-'^'*^ erg s^' AGNs, respectively). We therefore find that 
the average vLy{6ijm)/Lx ratios of Lx = 10'*^^'*^ erg s^' AGNs 
increase from ~ 3 at z ~ to 10 and ~ 43 at z = 0.5 — 1 and 
z = 1 — 2, respectively. By contrast, the average vLv (d/Jin) /Lx ra- 
tios of Lx = 10^3-44 erg s"' AGNs remain largely unchanged to 
z ~ 2 (i.e., vLv(6/jm)/Lx~ 1-2). Our data therefore suggest that 
deep IR observations can be used as a reliable method of identi- 
fying intrinsically bright (i.e., Lx= 104-'^44 gj.g g-l)^ potentially 
Compton-thick AGNs. However, we urge caution when using ei- 
ther the Lir/Lx or vLv(6/jm)/Lx ratios to locate less luminous ob- 
scured AGNs (i.e., Lx < lO^-' erg s^') unless other indicators of 
luminous AGN ac tivity are present (e.g., emission line spectra; see 
[Alexander et alj2008.) . 



Figure 6. Lm/Lx versus redshift for the stacked 70 (am undetected sources 
in the 'Restricted' sample of CDF-S X-ray AGNs, separated according to 
Lx (black squares). Vertical en'or bars indicate the range of Ljr/Lx pro- 
duced by assuming the various tracks in Figs. [1] and |4] Small, grey squares 
indicate the Ljr/Lx ratios of each X-ray AGN in the 'Restricted' sample 
derived from their 24 ^m flux densities, assuming the Average AGN and 
Average QSO SEDs for the Lx = lO^^-^a erg s"' andLx = IQ'^^-^* erg s"' 
AGNs, respectively (see ^3.3) . For those sources undetected at 24 \an we 
derive upper limits by assuming a 24 [om flux density limit of 50 [xTy. For 
comparison we include the individual and average Ljr/Lx ratios of the 
BAT/IRAS sample, separated according to Lx (vertical error bars indicate 
the standai'd deviation of Ljr/Lx of the BAT/IRS sample). Both the stack- 
ing analysis and the individual 24 [om detections clearly show a consistent 
increase in the Lir/Lx ratio of Lx = lO^^^^s ejg g-l AGN from z to 
z=l-2. 



age LiR of Lx = 104^-44 gj.g ^-i pj^^^ the z BAT/IRS 
sample is 1.2x10^^ Lq. Given the average X-ray luminosities of 
these AGNs (i.e., =3.8x10'*^ erg s"', 2.9x10"^^ erg s"' and 
3.2x 1043 erg s"' for the z «i 0, z = 0.5-1 and z = 1-2 AGNs, 
respectively), we constrain the change in Lir/Lx to a factor of 
~ I.8II0 between z = 0.5-1 and z = 1-2 and to a factor of 
~ l-2!'!og between z ~ and z = 1-2, both of which are consis- 
tent with no evolution; see Fig.|6] 

A number of recent studies have shown that the MIR emis- 
sion of AGN s can be used as a good proxy for their intrinsic lumi 



nosities (e.g. iKrabbe et"al] 1200 ll; iLutz et al.ir2004l ; iMaiolino et"al] 



I2OO7I; iHorst et alj|2008l) . Understanding how the X-ray-MIR re- 
lationship changes with redshift would be very useful for detect- 
ing potentially heavily obscured AGNs (e.g., iDaddi et al.l |2007| ; 
iFiore et al] |2008| ; iGeorgantopoulos et al. I I2OO8I) . However, previ- 
ous X-ray-MIR studies of AGNs have been limited to extremely 
high luminosity sources (i.e., Lx =3xl04^ erg s^^) at redshifts 
beyond z ~ 0.3. Using the results from our stacking analyses we 
can constrain the X-ray-MIR luminosity ratios (hereafter we fo- 
cus on vLv(6/jm)/Lx) of more typical, distant AGNs found in 
deep, multi-wavelength surveys. At z ~ the vLv(6/jm)/Lx ra- 
tios of the Average AGN and Average QSO SEDS are consistent 



^ We also note that the vLv(10.5^mi)/Lx ratios of the Average AGN and 
Average QSO at z~0 a re lai'gely consistent with the relationship presented 
in Krabbe et al] j200lh (i.e., vLv(10.5/;m)/Lx « 2.5 and 1.9, respec- 
tively). 



Table 1. The 30 X-ray detected sources with > 3a detection at 70 \xm. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


T06 


A03 


Name 


RA 


Dec 


z 


zflag 


log(Lx) (ergs s"') 


Nh (10^2 


^70 (AiJy) 


524 (AiJy) 


■^70 


LiR (i-o) 


Classification 


7 


312 


J033301. 8-275818 


53.2574 


-27.9718 


1.840 


0.6 


44.64 


3.57 


6370 ± 620 


640 ±11 


i\ ny 1 n r\t\ 

9.96 ± 0.99 


12.7 


A l_ M TOT A /^XT 

Abs. Non-BLACjN 


29 


240 


J033239. 0-275701 


53.1626 


-27.9505 


0.300 


0.9 


42.98 


5.32 


26600 ±2100 


o 'O o r\ 1 11 

3239 ± 11 


8.20 ± 0.64 


1 1.2 


A 1 M TOT A/^XT 

Abs. Non-BLACjN 


31 


230 


J033237. 9-275213 


53.1577 


-27.8704 


1.603 


3.0 


44.22 


1.79 


3920 ± 430 


777.4 ± 9.3 


5.05 ± 0.56 


12.3 


Abs. Non-BLAGN 


31 


230 


J033237. 9-275213 


53.1577 


-27.8704 


1.603 


3.0 


44.22 


1.79 


3920 ± 430 


o /I o o 1 o 

248.3 ± 9.3 


15.8 ± 1.8 


12.2 


A 1_ XT TOT A /^XT 

Abs. Non-BLACjN 


36 


204 


J033233. 1-274548 


53.1379 


-27.7634 


1.030 


0.5 


43.15 


1.73 


4020 ± 440 


477 ± 13 


8.42 ± 0.96 


1 1.9 


Abs. Non-BLAGN 


44 


173 


J033226. 6-274036 


53.1 108 


-27.6767 


1.031 


3.0 


43.89 


0.08 


8010 ± 770 


2107 ± 15 


3.80 ± 0.37 


12.2 


Unabs. BLAGN 


46 


163 


J033225. 3-274219 


53.1052 


-27.7054 


1.617 


3.0 


43.81 


1.08 


4580 ± 500 


485 ± 13 


9.4 ± 1.1 


12.6 


Abs. BLAGN 


51 


118 


J033217. 3-275221 


53.0719 


-27.8728 


1.097 


3.0 


44.01 


22.42 


A ^ Ar\ 1 ^ A f\ 

4340 ± 740 


259.2 ± 8.6 


16.7 ± 2.9 


1 1.7 


A 1_ XT TOT A /^XT 

Abs. Non-BLAGN 


52 


1 17 


J0332 17.2-274304 


53.0718 


-27.7178 


0.569 


3.0 


42.91 


0.04 


30500 ± 2300 


i ^ A^ 1 i r\ 

1247 ± 10 


^ A A 1 1 r\ 

24.4 ±1.9 


1 1.7 


TT L TOT A/^XT 

Unabs. BLAGN 


56 


88 


J033213. 3-274241 


53.0554 


-27.71 15 


0.605 


3.0 


43.31 


1.62 


r\/^or\ 1 o ^ f\ 

9620 ± 860 


665 ± 16 


14.5 ± 1.3 


11.3 


A 1_ XT T^T A /^XT 

Abs. Non-BLAGN 


69 


38 


J033201. 5-274138 


53.0063 


-27.6941 


0.850 


0.4 


43.40 


3.20 


1 /l C\r\ 1 f^^f\ 

1600 ± 620 


o c f\ 1 o r\ 

356.9 ± 8.9 


4.5 ± 1.8 


1 1.2 


A 1_ XT T^T A /^XT 

Abs. Non-BLAGN 


72 


30 


J033 158.3-275043 


52.9931 


-27.8453 


1.990 


0.5 


44.26 


7.77 


'^C\f\f\ 1 Af\f\ 

3020 ± 400 


^ f^y lO 1 r\ 1 

796.2 ±9.1 


3.79 ± 0.50 


12.4 


A 1_ XT TOT A /^XT 

Abs. Non-BLAGN 


77 


311 


J033301. 7-274543 


53.2570 


-27.7620 


0.622 


3.0 


42.53 


0.44 


5190 ± 540 


385.7 ± 9.1 


13.5 ± 1.4 


11.1 


Unabs. BLAGN 


78 


193 


J033230. 1-274524 


53.1256 


-27.7568 


0.960 


3.0 


43.06 


<0.19 


6960 ± 640 


AO A 1 1 O 

484 ± 12 


14.4 ±1.4 


1 1.7 


TT L T^T A/^XT 

Unabs. BLAGN 


98 


260 


J033244.3-275142 


53.1846 


-27.8618 


0.279 


3.0 


41.17 


<0.02 


91600 ± 6700 


7020 ±11 


1 O A C 1 f\ c\c 

13.05 ± 0.95 


11.5 


Unabs. Starburst 


122 


295 


J033257. 7-274549 


53.2404 


-27.7636 


2.100 


0.5 


43.47 


2.58 


4030 ± 460 








Abs. Non-BLAGN 


152 


302 


J033259.4-274859 


53.2476 


-27.8165 


1.280 


0.6 


43.80 


19.41 


3620 ± 430 


769 ± 1 1 


4.71 ± 0.56 


12.0 


Abs. Non-BLAGN 


155 


60 


J033208. 0-274240 


53.0335 


-27.71 1 1 


0.545 


3.0 


42.27 


3.59 


o o o r\ 1 A '\ f\ 

2880 ± 410 


455 ± 10 


o o 1 r\ r\o 

6.33 ± 0.92 


10.9 


A 1_ XT TOT A /^XT 

Abs. Non-BLAGN 


175 


284 


J033252. 0-274437 


53.2165 


-27.7438 


0.522 


3.0 


41.41 


<5.26 


51 IV ±580 


1213 ±11 


A ^ y 1 r\ A o 

4.76 ± 0.48 


1 1.2 


Abs. Starburst 


206 


109 


J033216. 3-273930 


53.0678 


-27.6586 


1.324 


3.0 


44.22 


0.13 


c 1 f\r\ 1 c A f\ 

5100 ± 540 


1376 ± 12 


o 1 1 r\ o r\ 

3.71 ± 0.39 


12.2 


TT L TOT A/^XT 

Unabs. BLAGN 


242 


285 


J03325 1 .9-274229 


53.2164 


-27.7083 


1.027 


3.0 


42.49 


0.72 


8150 ± 750 


500 ±11 


16.3 ± 1.5 


1 1.9 


T T L XT TO T A XT 

Unabs. Non-BLAGN 


253 


131 


J033220. 1 -274448 


CO r\o o 

53.0839 


-27.7467 


1.890 


1.9 


44.02 


73.51 


o 1 /^r\ 1 o or\ 

3160 ± 380 


327.20 ± 0.20 


9.6 ± 1.2 


12.4 


A 1_ XT T^T A /^XT 

Abs. Non-BLAGN 


268 


278 


J033249. 3-274050 


53.2056 


-27.6806 


1.222 


3.0 


44.10 


80.44 


A CiCf\ 1 C^f\ 

4850 ± 520 


O 1 O 1 1 

2127 ± 16 


o o o 1 r\ o c 

2.28 ± 0.25 


12.1 


A 1_ XT T^T A /^XT 

Abs. Non-BLAGN 


525 


129 


J033219. 9-274123 


53.0828 


-27.6899 


0.229 


3.0 


41.09 


0.00 


14000 ± 1200 


2932 ± 17 


A 1 r\ A t 

4.79 ± 0.41 


10.7 


T T L XT TO T A XT 

Unabs. Non-BLAGN 


538 


65 


J033208.6-274649 


53.0359 


-27.7803 


0.310 


3.0 


41.23 


0.51 


12900 ± 1100 


973 ± 14 


13.2 ± 1.1 


10.7 


Unabs. Non-BLAGN 


553 


292 


J033256.8-275318 


53.2366 


-27.8885 


0.366 


3.0 


41.27 


<0.17 


8520 ± 770 


612.3 ±7.1 


13.9 ± 1.3 


10.7 


Unabs. Starburst 


567 


236 


J033238.9-274732 


53.1620 


-27.7925 


0.460 


3.0 


41.05 


<0.16 


6370 ± 600 


730.9 ± 8.9 


8.72 ±0.83 


11.1 


Unabs. Starburst 


575 


115 


J033217.2-274922 


53.0715 


-27.8230 


0.340 


3.0 


41.02 


1.90 


3290 ± 370 


221.5 ±9.5 


14.8 ± 1.8 


10.2 


Abs. Non-BLAGN 


577 


224 


J033236.3-274933 


53.1512 


-27.8258 


0.547 


3.0 


41.42 


<0.14 


10470 ± 890 


664 ± 10 


15.8 ± 1.4 


11.2 


Unabs. Starburst 


608 


318 


J033304.0-275027 


53.2665 


-27.8409 


0.890 


3.0 


43.86 


150.00 


3090 ± 410 


578.0 ± 9.3 


5.34 ±0.71 


11.6 


Abs. Non-BLAGN 


646 


266 


J033245.2-274724 


53.1884 


-27.7903 


0.438 


3.0 


41.10 


<0.19 


11700 ± 1600 


507.2 ± 8.6 


23.1 ± 3.2 


11.0 


Unabs. Starburst 



Notes: (1) Source index from lTozzi et al.l f2006h . (2) Source index from lAlexander et al.l i2O03l) . (3) Source name. (4-5) X-ray position (J2000). (6) Source redshift. (7) Redshift flag. (8) Intrinsic (i.e., corrected for 
absorption) 2 — WkeV luminos ity. (9) Colu mn density in units of lO^^cm^^. (10) 70 \Mn flux. (1 1) 24 |-im fl ux. (14) 70/24 |jm flux ratio. (13) Ljr derived from the 70 [om flux using the method described in the text. 
(14) Classification from , Bauer et 1l H2004h . Information in cols. (l),(3)-(9) is taken from lTozzi et alj j2006l) . 



Table 2. Descriptions and average properties of the AGN and SB stacks. 











(3) 






(4) 






(5) 






(6) 






(7) 






(8) 




(1) 


(2) 






N 






Mean z 




log(Lx) (ergs s 






Stacked S70 (pJy) 




Stacked ^24 0"Jy) 




LlR (L0) 




Index 


Description 




All 


Det. 


Und. 


All 


Det. 


Und. 


All 


Det. 


Und. 


All 


Det. 


Und. 


All 


Det. 


Und. 


All 


Det. 


Und. 


1 






14 




3 


0.46 


0.44 


0.47 


41.41 


41.26 


41.49 


9960 =b son 


22400 ± 1200 


630 ± 110 


736.8 


1592.9 


94.7 


11.0 


11.3 


10.1 


2 


z < 1 = 10^ 


1-42 


25 


3 


22 


0.43 


0.29 


0.45 


41.60 


41.12 


41.64 


1540 ± 160 


10070 ± 550 


380 zt 170 


235.7 


1183.5 


106.4 


10.4 


10.8 


<10.5 


3 


z > 3 




12 





12 


3.52 




3.52 


44.31 


- 


44.31 


280 ±180 




280 ±180 


127.5 


- 


127.5 








4 


BLAGN 




34 


6 


28 


1.54 


1.02 


1.65 


43.94 


43.74 


43.97 


2050 ± 150 


10050 ± 450 


330 ± 160 


419.3 


957.3 


304.1 


12.0 


12.3 


<11.9 


5 


z = 0.5 - 1 


= l()«-44 < 1022 


32 


3 


29 


0.74 


0.72 


0.74 


43.10 


42.88 


43.12 


2180 it 170 


14200 ± 830 


940 ± 170 


301.6 


820.2 


247.9 


11.2 


11.7 


10.9 


6 


7 = 1 - 2Lv = 


JQ42-44 ^ jq22 


19 


2 


17 


1.38 


1.03 


1.42 


43.39 


43.61 


43.35 


1230 ± 160 


8080 ± 540 


420 ± 170 


265.0 


1156.6 


160.1 


11.7 


12.2 


<11.7 


7 


z = 0.5 - 1 


= 10«-44 ^ 1022-23 


31 


3 


28 


0.75 


0.67 


0.76 


42.99 


43.20 


42.96 


750 ± 140 


4700 ± 380 


330 ± 150 


180.0 


494.0 


146.3 


10.8 


11.4 


<10.9 


8 


z = 1 - 2 = 


10«-*'' Nh = 10^2-23 


32 


2 


30 


1.38 


1.32 


1.39 


43.35 


43.59 


43.33 


680 ± 140 


4300 ± 330 


440 ± 150 


174.4 


480.3 


154.0 


11.4 


12.3 


<11.7 


9 


z = 0.5 - 1 La- 


= 1042-44 < io23 


63 


6 


57 


0.74 


0.69 


0.75 


43.05 


43.07 


43.05 


1440 ± 110 


9450 ± 460 


600 ± 110 


241.7 


657.1 


198.0 


11.0 


11.5 


10.7 


10 


7 = 1 — 2 Ly = 


JQ42-44 ^ jo23 


51 


4 


47 


1.38 


1.18 


1.40 


43.36 


43.60 


43.34 


920 ± 110 


6190 ± 320 


470 ± 120 


208.1 


818.5 


156.2 


11.5 


12.3 


11.2 


11 


z = 0.5 - 1 Lx 


= 10«-« Af/, < 10^3 


44 


3 


41 


0.74 


0.58 


0.75 


42.63 


42.65 


42.63 


1390 ± 130 


12850 ± 820 


550 ± 130 


167.8 


688.3 


129.7 


11.0 


11.4 


10.7 


12 


z=l-2Lx = 


10«-43 Njj < 1023 


19 


1 


18 


1.35 


1.03 


1.37 


42.68 


42.49 


42.69 


1110 ± 170 


8150 ± 750 


720 ± 170 


153.0 


555.6 


130.6 


11.7 


11.9 


11.4 


13 


z = 0.5 - 1 Lx 


= 10«-4* < 1023 


19 


3 


16 


0.77 


0.81 


0.76 


43.43 


43.28 


43.46 


1510 ± 200 


6060 ± 410 


660 ± 230 


413.1 


625.9 


373.2 


11.1 


11.8 


<11.2 


14 


z=l-2Lx = 


lO^-"* Nh < 1023 


32 


3 


29 


1.40 


1.23 


1.42 


43.53 


43.72 


43.51 


870 ± 140 


5540 ± 340 


390 ± 150 


240.9 


906.1 


172.1 


11.5 


12.3 


<11.7 


15 


SiA = SO/jJy — 


150;j7y 


68 





68 


1.40 




1.40 


43.73 




43.73 


260 ± 90 




260 ± 90 


113.2 




113.2 


<11.4 




<11.4 


16 


S24 = ISQ/jJy- 


- 500ijJy 


68 


10 


58 


1.37 


1.00 


1.44 


43.75 


43.52 


43.78 


1263 ± 97 


4420 ± 180 


720 ± 110 


285.9 


456.9 


256.4 


11.7 


11.9 


11.5 


17 


524 ^ SOOfjJy 




20 


13 


7 


1.10 


1.01 


1.27 


44.02 


44.01 


44.04 


7040 ± 230 


10120 ± 300 


1330 ± 320 


1113.4 


1199.7 


953.0 


12.2 


12.3 


11.6 


18 


S24 ^ 50^(/>' 




156 


23 


133 


1.35 


1.01 


1.41 


43.79 


43.86 


43.77 


1570 ± 66 


7640 ± 190 


520 ± 70 


316.7 


876.8 


219.9 


11.7 


12.1 


11.3 



Notes: (1) Stack index. (2) Stack description. (3) The numbers of CDF-S AGNs/SBs in each stack. (4) Mean redshift of the sources in each stack. (5) Logarithm of the mean, absorption corrected X-ray luminosity 
of the sources in each stack. (7) Stacked 70 (jm flux density per source. (8) Stacked 24 pm flux density per source. (9) Infrared luminosity per source (over the 8-1000 [jm interval), derived from the 70 |jm flux 
density. Upper limits include the uncertainty on the S70-L1R correction. We do not calculate Lir for stack (4), z > 3 AGN, as the mean redshift is beyond the range of the tracks derived from the BAT/IRS sample. 



3.3 24 Properties 

The results on the Ljr/Lx ratios in i]3.2l were based on stacking 
analyses, which could be biased by a few bright 70 \im undetected 
sources. Since the majority of the CDF-S X-ray AGNs are indi- 
vidually detected at 24 [xm, we can perform a complementary test 
of our results that does not rely significantly on stacking analyses. 
The advantage of this approach is that we can assess the range of 
LiR for the X-ray AGNs, although there can be considerable uncer- 
tainty in the conversion between 24 [im flux density and Lir for 
individual objects; see Fig.jS] However, as we show in Fig.jS] pro- 
vided the average SED is not predominantly starburst dominated, 
then we can predict the average Ljr/Lx ratio, on the basis of the 24 
^.m data, to within a factor of ~ 3; on the basis of Fig. ID the aver- 
age SEDs are AGN dominated, and would be increasingly so if the 
stacked 70 \un flux is dominated by bright undetected sources. Fur- 
thermore, so long as the average SEDs of the z = 0.5-1 AGNs are 
the same as the z = 1-2 AGNs, then we can accurately assess the 
relative change in Lir/Lx between these redshifts; we provide evi- 
dence in §4. 1 that this does appear to be the case. To convert from 
24 \im flux densities to Ljr we use the same SEDs as assumed 
when calculating Lir from the stacked 70 [im flux densities (i.e.. 
Average AGN and Average QSO for the Lx= lO^^^^^* erg s^' and 
jq43-44 gj.g g-l ^Qjvjj^ respectively; see ^3.2.2\ . For those AGNs 
that are undetected at 24 ^im we calculated upper limits on Ljr as- 
suming S24 = 50 (i-Jy; however, we note that this assumption has 
a small effect on our average results (i.e., Lir/Lx only changes by 
~ 35% even if we assume extreme upper limits of ^24 = (i-Jy). 
The average Lir/Lx ratios derived from the 24 |j.m fluxes are given 
in Table 1. 

We include in Fig. [6] the individual Ljr/Lx ratios calculated 
using the Lir values derived from 24 \im fluxes of the 70 |i.m 
undetected AGNs in the 'Restricted' sample. For both the Lx= 
jQ42-43gj.g 10'*^"'*'*erg s^i AGNs the mean Lir/Lx 

ratios increase with redshift by factors of Ri4.5 and ~1.8, respec- 
tively. These results are in qualitative agreement with those ob- 
tained from the 70 [i-m stacks and indicates that the stacks are not 
dominated by a few extreme sources; these results are also con- 
firmed by the median Lir/Lx ratios presented in Table [3] How- 
ever, we note that there is a factor of ~2 offset between the Lir/Lx 
ratios derived from the 70 \im. stacks and individual 24 \im. flux 
densities. This is caused by the difference between the S70/524 
flux ratios of the 70 |j.m undetected stacks and the assumed SEDs 
(see Fig. ^ and disappears if we assume the SEDs of the two 
individual BAT/IRS AGNs that most closely reproduce (at both 
redshift bins) the Sh)/S2A flux ratios of the 70 \im undetected 
Lx= 1042-43 gj.g ^o43-44 gj.g ^-1 ^j^^j,^. ^qq ^275 and 

Mrk 290, respectively. 



4 DISCUSSION 

We have provided multiple lines of evidence for an increase in the 
Lir/Lx ratio for AGNs over the redshift range z = — 2, with Lx= 
jq42-43 gj-g g-l AQNs ~5 times more infrared luminous at z = 1 — 
2 than at z = 0.5 — 1 (and ~ 20 times more infrared luminous than at 
z ~ 0). The evidence for an increase in Lir for Lx= 104-'-44 gj.g g- 1 
AGNs is less conclusive: a factor of ~2 between z = 0.5 — 1 and 
z = 1 — 2 and no change from z ~ and z = 1 — 2. 

These results provide new insight into the production of IR 
emission from distant AGNs and lead to a number of questions, 
which we address below. 



4.1 What is driving the increase in Lir? 

The Lir/Lx ratio for the Lx= lO^^-^S gj-g ^-i ^.GNs at z = 1 - 2 
is higher than that found for lower redshift X-ray AGNs; see Table 
[3l However, the Lir/Lx ratio for the Lx= lO^^^^S gj-g g-l AGNs 
is broadly consistent with z ~ 2 submillimeter emitting galax- 
ies (S MGs) hosting AGN activity (Lftp /Ly~250; [Alexander et all 
I2OO5I) . On the basis of sensitive Spitzer-\KS spectroscopy, the 
large Lfir/Lx ratio from AGN-hosting SMGs appears to be due 
to intense star-formation activity, with an average contribution 
to LpiR from AGN activity of a:il O% (e.g.. Pope et al. 200^; 
iMenendez-Delmestre et al .I200ll200^ . Can the increase in the av- 
erage Lir/Lx ratio of the z = 1 — 2 AGNs also be due to increased 
star-formation activity? 

There is no clear relative offset between the stacked 
S70/524 flux ratio and the AGN-dominated track for the Lx= 
]^q42-43 gj-g g-1 AGNs over the redshift range z = 0.5 — 2, ap- 
pearing to suggest that the relative AGN/SB contribution to Lir 
has not changed with redshift. If we derive the relative AGN/SB 
contribution to the stacked S-jq/Sia using the average AGN and 
SB tracks, then we estimate ai60% and fa70% of the AGNs with 
Lx= 1042-43 gj-g g-i have AGN-dominated SEDs at z = 0.5 - 1 
and z = 1—2, respectively; by comparison, the stacked 570/524 
of the Lx= 1043-44 erg s"' AGNs suggest f»100% have AGN- 
dominated SEDs over the full redshift range. However, since it is 
difficult to unambiguously determine the relative AGN/SB contri- 
bution to Lir from the S-1Q/S2A flux ratios at z > 1.5 (see Fig.|4]l, 
this result should be considered tentative^ Therefore, we explore 
below the implication of our results assuming both an increase in 
star formation and AGN activity. 

If the observed increase in Lir is attributed to star-formation 
it would imply a significantly higher (i.e., a factor of ~5) ratio be- 
tween star-formation and black hole growth at z = 1 — 2 compared 
to z = 0.5 — 1, and a factor of ~ 13 increase between z ~ and 
z = 1—2. Whether this has significant implications for the black- 
hole-bulge mass relationship depends on the location of the star- 
formation in the host galaxy and on the AGN fraction (i.e., the 
duty cycle of black-hole growth). For example, using Eqn. 4 of 
lKennicuttl ( |l998,) to derive the star-formation rate from Lir and de- 
riving the mass accretion rate from Lx, we estimate that the ra- 
tio between the average star-formation and mass accretion rate is 
~500 in the z = 1 — 2 AGNs, which is consistent with the black- 
hole to bulge mass ratio observed in the local uni verse (i.e., ~800; 
iMcLure & Dunlotj|2002l ; lMarconi & Huntll2003i) . This would ap- 
pear to suggest a closer relationship between black-hole and stellar 
growth at z = 1—2 than found at z = 0.5 — 1, where the average 
star-formation and mass accretion rate would be ~100 (and ~30 
at z ~ 0), based on the same assumptions as above. However, these 
different results could be reconciled if the AGN fraction is higher at 
higher redshifts or if the majority of the star formation at higher red- 
shifts is occurring in the galaxy disk rather than the galaxy bulge, 
or vice-versa. 

If, however, the increase in the Lir/Lx ratio is due to the 
AGN, it would imply that a larger fraction of the intrinsic emis- 



We note that we also find a similar fraction 30%) of Lx= 
1042-« erg AGNs at z = 1 - 2 have IRAC 3.6-8.0 Mm colours sug- 
gesting AGN-dominated SEDs (based on IStem et aLlboOSl) as that found 
for Lx= 104^-4^ erg s-' AGNs at z = 0.5 — 1. This also suggests that the 
relative AGN/SB contribution forLx= iC^-^a grg s-' AGNs is constant 
with redshift but we caution that this result is based on shorter-wavelength 
data than used to derive Ljr. 



Table 3. Average properties of the CDF-S X-ray AGNs in tlie 70 \im undetected, 'Restricted' sample (see i|2.1.U . Naii is the number of AGNs detected in each 
bin and Nj/,, is the number of these that are detected at 24 |am. z is the mean redshift of the AGNs in each subsample and log(Lx) and log(Lx) indicate their 
mean and median Lx- We provide the mean Ljr/Lx ratios derived from the 70 [am and 24 |jm stacks (assuming the Average AGN and Average QSO SEDs 
for the Lx= 10'*-"'*' erg s"' and Lx= 10'*^^'*'* erg s"' subsamples, respectively) as well as the mean and median Ljr/Lx ratios derived from the individual 
24 [am flux densities. Finally, we provide the mean and median Lir/Lx derived from the individual 24 \mi flux densities assuming the two BAT/IRAS tracks 
that most closely match the Sio/S^a flux ratios of the Lx= 10'*^^'*^ erg s^' and Lx= 10'*^"'*'' erg s"' subsamples at both z = 0.5 — 1 and z = 1 — 2: NGC 1275 
and Mrk 290, respectively. 
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sion from the accretion disk is reprocessed by dust. This increase 
in the Lir/Lx ratio may therefore be due to larger AGN dust cover- 
ing factors at higher redshifts; there is tentative evidence that this is 
indeed the case from measurements of the o bscured to unobscured 
ratios of AGN s in deep X-ray surve vs (e.g. iLa Franca et al.ll2005l : 
lHasingej2008ll . Results presented in lHasingeil ( l2008h predict a fac- 



tor of ~2 increase in the ratio of obscured-to-unobscured AGNs 
between z = and z = 1 — 2, significantly lower than what we mea- 
sure here. However, this is based on a sample of AGNs covering a 
broader range of Lx than focused on here and there is evidence that 
any increase in the dust covering fraction will be weaker in more 
X-ray luminous AG Ns (e.g. Ueda et al. 2003; Akvlas et al. 200i; 
iTreister et al J2008l ; and supported by the reduced Ljr/Lx ratios of 
the Lx= 10^^^^^ erg s^' AGNs reported here). Therefore, on the 
basis of the X-ray survey results, it is plausible that at least some of 
the increase in Ljr/Lx ratio is due to a larger dust covering factor at 
z ~ 1 — 2 than seen at lower redshifts. An increase in dust-covering 
factor with redshift is predicted by a number of theoretical mod- 
els, which suggest that AGNs undergo early growth during a hid- 
den phase before expelling their obscuri ng gas and dust, revealing 
a luminous, tmobscured q uasar ^Silk & Reeslll99^ ; ISpringel et al.l 
12001 iHopkins et alfeOOQ) . 

More direct constraints on the origin of the increase in the 
Ljr/Lx ratio will be placed using the Herschel Space Observatory 
(see 343} and Spitzer-IKS, for the fraction of X-ray AGNs that have 
Spitzer-IKS data (J. R. Mullaney, in preparation). 



4.2 What is the contribution of AGNs to the cosmic IR 
background? 

We can use the results of our stacking analysis to place constraints 
on the AGN contribution to the 70 \im. background. If we stack all 
of the 251 X-ray AGNs, irrespective of whether they are individ- 
ually detected at 70 \im we obtain an average 70 \xm flux density 
of 1040±80 (i-Jy, and therefore an integrated 70 vim flux density of 
260±20 mJy, over the 391.3 arcmin^ of the CDF-S field. On the 
basis of the analysis of the I R background in both GOODS fields 
and the Lockman Hole (e.g.. lDole et alJuOOfil) , this corresponds to 
>5% of the average resolved 70 \\m background. This constraint 
is a lower limit since (1) the X-ray observations will not have iden- 
tified the most heavily obscured luminous AGNs in this field, (2) 
the 70 \im field is too small to include X-ray bright AGNs, and (3) 
the X-ray observations are not sensitive to the lowest-luminosity 
AGNs. 

On the basis of iTozzi et al. I ( l2006l) , we would predict ~80% 
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Figure 7. Top: Predicted fractions of CDF-S X-ray AGNs detected in the 
upcoming HGOODS deep infrared survey to be undertaken by the Herschel 
Space Obsen'atory . The expected observed 100 |jm flux densities are cal- 
culated by extrapolating the Spitzer-MWS 24 ytm flux densities along the 
Average SB-dominated and Average AGN-dominated SEDs, derived from 
the BAT/IRS sample of AGNs. If all the X-ray AGNs were to have SB- 
dominated IR SEDs, we would expect to detect at least as high a fraction 
at 100 \an as we currently detect at 24 \mi. Also shown for comparison are 
the fractions of CDF-S AGNs detected by MIPS at 70 |am in each redshift 
bin Bottom: The number of CDF-S X-ray AGNs in each redshift bin. 



of Lx> lO'^'erg s^^, N^< 10^3 AGNs in the CDF-S (i.e., to z ^5) 
to be X-ray undetected. If we assume that they have the same IR 
SED as the X-ray detected AGNs then they increase our estimate 
of the resolved 70 \xm backgrou nd by a factor of five to >25%; 
from a study of nearby sources, iLutz et al. 

I ( l2004l) finds that the 
IR emission of AGNs is not significantly depressed in the most 
heavily obscured sources (we confirm that this is also the case for 
the BAT/IRS sample). Since X-ray surveys are insensitive to the 
most heavily obscured AGNs (i.e., Compton-thick sources with 
A'h> 10^^ cm^^), the overall contribution to the 70 yt-m back- 
ground could potentia lly double to > 10% (e.g., lDaddi et al.f2007l ; 
[Alexander etalj2008l) . 

To estimate the contribution to the 70 ^.m background from 
bright AGNs that are too rare to lie in the small CDF-S field 
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Figure 8. AGN/SB diagnostic powers of the PACS filters on board Her- 
schel Space Observatory. The shaded regions indicate the range of expected 
flux ratios of the AGNs in the BAT/IRS sample at the redshifts shown. Solid 
lines indicate the expected flux ratios of the BAT/IRS Average AGN and 
Average SB. 



would require the same analyses as performed here to be made 
on a shallower survey that covers a larger area of the sky. As it 
stands, the most accu rate 70 [im analysis in a larger field is that of 
IPapovich et al] j2007t) of the E-CDFS. They find that bright AGNs 
(i.e., Lx> lO'^'^erg s"') contribute only 7x10"'* MJy sr"' to the 
IR background at 70 (xm, corresponding to an insignificant fraction 
(<^C 1%) of the total IR background at this wavelength. However, 
the E-CDFS field is not a great deal larger than the CDF-S field 
and a full study of the contribution of X-ray bright AGNs to the IR 
background will require the analysis of the much larger COSMO^ 



and/or Booteo fields at 70 [im. 



4.3 What advances may we expect for deep surveys with 
Herschel"! 

Our study of the FIR properties of X-ray detected AGNs allows 
us to place constraints on the properties of AGNs detected in deep 
Herschel fields at > 70 \im. We anticipate that the improved sen- 
sitivity of the Herschel-PACS instrument over Spitzer-MIPS will 
result in a significant increase in the fractions of X-ray AGNs that 
will be detected at FIR wavelengths, which will provide more direct 



' URL: |http://cosmos.astro.caltech.edu/| 
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insight into the processes driving the increase in Lir at z = 1 — 2. 
Depending on the assumed SED we predict that approximately 45- 
60%, 50-70%, 50-65% and 45-65% of the CDF-S X-ray sources at 
z<I,z=l — 2,2 — 3 and 3 — 4, respectively, will be detected at 
100 |J.m in the proposed ultradeep HGOODS survey (limiting flux 
= 0.6 mjy; see Fig. [7}. These estimated detection levels are calcu- 
lated by extrapolating the 24 ii.m flux densities along the Average 
SB-dominated and Average AGN-dominated SEDs (derived from 
the BAT/IRS sample). 

In Fig. [8] we present the results of passing the IR spectra of the 
BAT/IRS sample through the various PACS filter response curves. 
Our analysis shows that, although the S100/S24 flux ratio will not 
push our constraints on the relative contributions from AGN activ- 
ity and star formation to higher redshifts (due to the influence of 
the MIR spectral features that shift into the observed 24 [im band 
at high redshifts), the Sioo/Sjo and Si^q/Sjo flux ratios will "take 
over" and provide this information out to a redshift of ~ 5 (i.e., the 
redshift of the most distant AGNs currently detected in the deep- 
est X-ray surveys). Deep, infrared observations undertaken by Her- 
schel will therefore provide a method of identifying a significant 
proportion of the long sought after population of Compton-thick 
AGNs out to z ~6. 



5 SUMMARY 

We have investigated the MIR and FIR properties of X-ray de- 
tected AGNs and SBs in the CDF-S using deep 24 |j.m and 70 [im 
observations undertaken by the Spitzer Space Telescope as part 
of the GOODS and FIDEL legacy surveys. Out of the 266 X- 
ray AGNs/SBs, 30 (sill%) and 172 (^65%) have counterparts at 
70 vim and 24 [im respectively, with a bias towards the detection 
of SBs and BLAGNs at 70 |i.m. As the majority of CDF-S X-ray 
AGNs are undetected at 70 ^.m we rely on stacking analysis to 
measure their average MIR and FIR properties. We compare the 
IR properties of the CDF-S AGNs/SBs with those of a sample of 
local AGNs that have X-ray properties (i.e., Lx and A'h) covering 
the same range as the CDF-S AGNs (i.e., the BAT/IRS sample). 
In the following points we summarise the main conclusions of this 
study: 



(i) We find strong evidence at both 24 \xm and 70 \im that 
the average infrared luminosity of Lx= 10'*^"'*^ erg s"' AGNs 
at z = I — 2 is significantly higher than those at z = 0.5 — I and 
z ~ (i.e., by a factor of ~ 5 and ~ 20, respectively). This dif- 
ference corresponds to a factor of 4.7^2*0^ ^^'^ ^'^■'^-I'e increase 
in Lir/Lx, respectively. This large increase in Ljr is not, how- 
ever, seen for AGNs with higher X-ray luminosities (i.e., Lx= 



10-* 



erg s ). We therefore argue that deep IR observations 



can be used as a reliable method to identify intrinsically bright 



(i.e.., Lx= 10' 



,43-44 



erg s ) Compton-thick AGNs, but urge cau- 



tion when using this ratio to locate less luminous obscured AGNs 
(i.e., Lx < lO"*"* erg s"'), unless other indicators of AGN activity 
are present. 

(ii) Due to the low numbers of AGNs detected at 70 [im we are 
unable to establish what process is driving this increase in the av- 
erage Lir; however, both increased star-formation and/or increased 
AGN dust covering factors are likely candidates. If the former, then 
the z = 1 — 2 epoch may represent a period of rapid growth of the 
bulge to black hole mass ratio. However, there is tentative evidence 
from X-ray observations that the dust covering fraction is, indeed, 
higher at large redshifts. We predict that forthcoming deep surveys 
to be carried out by the Herschel Space Observatory will enable us 



to resolve what process is driving the increase in Lir. See sections 
[3:2:211331141] 

(iii) On average, more X-ray luminous CDF-S AGNs have lower 
S-1Q/S2A flux ratios. Based on the infrared properties of a sample of 
local AGNs with similar X-ray properties (i.e., Lx and A'h), we con- 
clude that more X-ray luminous CDF-S AGNs have IR SEDs that 
are more AGN-dominated (rather than SB dominated). See section 

[mi 

(iv) Despite measuring a large increase in Lir among z = 1 — 2 
AGNs (see point[(i)]above), we find that the X-ray detected AGNs 
in the CDF-S contribute only ~5% of the 70 vim background. How- 
ever, if we extrapolate this fraction to take into account those ~80% 
of AGNs that are undetected in X-rays, we estimate that ~25% of 
the 70 vim background is attributable to AGNs. See section[4j2l 

(v) We anticipate that the undertaking of deep, FIR surveys by 
the Herschel Space Observatory will allow us to detect between 
40% and 75% of the X-ray detected AGNs in the 1 Ms CDF-S, 
depending on whether the infrared SEDs are predominantly AGN 
or SB-dominated, respectively. Furthermore, the FIR diagnostics 
used in our analysis can be directly applied to the R4CS filters on 
board Herschel to enable us to discriminate between AGNs with 
SB-dominated or AGN-dominated infrared SEDs. See section [431 
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APPENDIX A: THE Swift-BAT COMPARISON SAMPLE 

To aid in the interpretation of the Spitzer MIPS data used to char- 
acterise the CDR-S X-ray sample, we selected a local sample of 
AGNs from the Swift-BAT survey with Spitzer-IRS low-resolution 
spectroscopy (5.2-38 vim). The Spitzer-IRS data allows us to ac- 
curately distinguish between AGNs with AGN-dominated and SB- 
dominated IR SEDs on the basis of the MIR spectral features in 
the Swift-BAT AGNs. The Swift-BAT survey provides the ideal lo- 
cal AGN comparison sample since the sensitivity of the BAT tele- 
scope to hard X-ray photons (14-195 keV) provides an almost 
absorption-independent selection of AGNs (to A'h~ 10^^ cm^^). 
Indeed, the range of X-ray luminosities and absorbing column den- 
sities of the Swift-BAT AGNs are comparable to those of the CDR-S 
AGNs (see Tueller et al. 2008 and Winter et al. 2009); see §2.2.1. 

Thirty-six of the 104 Swift-BAT AGNs with g ood-quality 
I^2-m kfV snd N}i constraints from i Winter et alj ( [2009 ) or 
iBassani et all ( 1 19991) have both Short-Low (SL; 3.6" x 136.0"; 
R ~ 60-127) and Long-Low (LL; 10.5" x 360"; R ^ 57-126) 
Spitzer-TRS spectroscopy, providing full coverage at 5.2-38 ^m. 
All 36 BAT/IRS AGNs were observed in staring-mode with two 
nod positions, which are required for background subtraction. Ba- 
sic Calibrated Data (BCD) images were combined and cleaned as 
described in Goulding & Alexander (2009). Differing nod positions 
were subtracted from each other to produce background subtracted 
images before extracting the spectra using SPIC^lB The objects 
were sufficiently bright and the observations were sufficiently short 
that there was no significant effect from latent charge build-up on 
the detector. Long-slit orders were clipped using the wavelength 
trim ranges given in Table 5.1 of the Spitzer-IRS Observers man- 
uaQ The continuum of each source in every order was fitted using 
a first or second-order polynomial. The spectra were then matched 
to give a single continuous spectmm for each source. Rlux calibra- 
tion was carried out using the latest available Spitzer-IRS calibra- 
tion files (version 17.2) and is largely consistent with archival IRAS 
25 [im flux densities (to within a factor of 2.5 in ~97% of cases). 
The BAT/IRS AGNs have absorption corrected 2-10 keV luminosi- 
ties spanning 4.3x 10^' erg s^' to 1.9x 10^^ erg s^' and absorbing 
column densities in the range 2.5 x lO^O cm-2 to 2.2x10^4 cm 

" URL:lhttp://ssc. spitzer.caltech.edu/postbcd/spice. html" 

The Spitzer-IRS Observers Manual is available at 
|http://sscTspitzer.caltech.edu/irs/dh/| 
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Figure Al. IRS spectra of the 25 BAT/IRS AGNs that we classify as having AGN-dominated IR SEDs on the basis of the lack of a prominent PAH feature at 
6.2 |am. Also indicated ai'e the PAH features at 7.7 \mi and 1 1.4 \xm although these were not used to distinguish between AGN and SB dominated systems. 



(i.e., roughly the same as those C DF-S X-ray AGN /SBs with Lx 
and Nil measurements published in iTozzi et aljioO^) . 

To determine whether the S-iq/S2A flux ratio can distin- 
guish between SB and AGN dominated systems, we classified the 
BAT/IRS sample into those objects with a) a prominent PAH fea- 
ture at 6.2 \x.m, which we assume is a go od indicator of a str ong 
starburst component (SB-dominated; e.g.. lOenzel et alj|l998l : 11 
objects), and b) a power-law MIR to FIR continuum with no sign 
of any PAH feature at 6.2 \im, which we assume indicates an AGN- 
dominated object (25 objects). The 6.2 \x.m PAH feature is preferred 
over other lines as it lies in a region of the IR spectrum that is 
largely free of other strong spectral features, in particular the sili- 
cate absorption/emission band at 9.7 ii.m. We show the 5p;rzer-IRS 
spectra of the BAT/IRS sample in Fig. Al (AGN-dominated sys- 
tems) and Fig. A2 (SB-dominated systems). 

To calculate the expected 570/S24 flux ratio tracks each of 
for the 36 BAT/IRS AGNs as a function of redshift, we shift each 
spectrum by a factor of 1 + z in wavelength (where z is increased 
from 0.25 to 2.5 in steps of 0.025) and pass the resulting spec- 



trum through the MIPS 24 |j.m and 70 \im filters. To determine the 
observed 70 \im fluxes at z < 0.75 requires knowledge of the IR 
SED beyond the wavelength coverage of the Spitzer-[RS> spectra 
(i.e., > 38 |J.m). To provide this longer wavelength coverage, we 
extrapolated the longest wavelength bin of the IRS spectra to the 
60 [xm and 100 \im flux densities from the IRAS Point Source and 
Faint Source Catalogue^ Thirty-two of the 36 AGNs with IRS 
data have well constrained IRAS 60 ii.m and 100 \im. fluxes. The 
remaining 4 (all classed as AGN dominated) have IRAS 60 |i.m or 
100 (im flux densities that are flagged as upper limits in the IRAS 
catalogues. However, including these faint sources in our average 
AGN SED by assuming the upper limits as detections has no effect 
on any of our results. In Fig. I A3 1 we show the 570/524 flux ratios 
derived from the BAT/IRS sample as a function of redshift. There 
is a clear separation between AGN and SB dominated systems at 
z < 1.5, showing that the S^o/S24 flux ratio can be used as an effi- 



Available from:lhttp://irsa.ipac.caltech.edu/applications/Gator/index.html| 
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Figure A2. IRS spectra of the 1 1 BAT/IRS AGNs that we classify as having SB-dominated IR SEDs on the basis of the presence of a prominent PAH feature 
at 6.2 (jm. Also indicated are the PAH features at 7.7 yim and 1 1.4 (jm although these were not used to distinguish between AGN and SB dominated systems. 



cient method to determine whether AGN or star-formation activity 
dominates the IR emission of X-ray detected AGNs. 

Calculating Lir for high redshift sources using the 24 |j.m flux 
density alone is susceptible to systematic errors caused by spectral 
features that are redshifted into this waveband at z > 0.7. We can 
use the BAT/IRS sample to establish the uncertainties in deriving 
LiR from either 70 \im or 24 [xm data. In the lower panels of Fig. 
IA3l we show the expected observed frame 70 vim and 24 \im fluxes 
for the BAT/IRS sample, if observed at z=0.25-2.5 and normalised 
to LiR = 10^^ Lq (the normalisation factor, lO'^ Lq/Ljr, is derived 
by calculating Ljr from the IRAS flux densities usi ng the equations 
presented in Table 1 of Sanders & Mirabel|[T993) . Depending on 
the shape of the IR SED, sources of the same Lir can differ in their 
observed 24 |i.m flux density by up to a factor of ~12 at z > 1, 
compared to factor of only ~3 at 70 (im. This range of 70 \im flux 
density is almost independent of redshift out to z = 2.5, indicating 
that is up 4 times more accurate for determining Lir when no 
information on the SED shape is available. 
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Figure A3. Top.The expected 570/524 flux ratio of the AGNs in the BAT/IRS AGNs over the redshift range z = 0.25 — 2.5. Middle: The expected observed 
70 [am flux densities of the BAT/IRS AGNs over the redshift range z = 0.25 — 2.5, each normalised to Lir= lO'^ Lq. Bottom: The expected observed 24 yrni 
flux densities of the BAT/IRS AGNs over the redshift range z = 0.25 — 2.5, again normalised to Lir= lO'^ Lq. We show the ratios and flux densities as 
expected if the AGNs were to be observed atz = 0.25 — 2.5. Solid, red line s refer to AG N-dominated and dashed blue lines to SB-dominated AGNs. The green 
shaded areas indicate the range of ratios and flux densities spanned by the lBrandl et alj ^200© starburst galaxies over the redshift range z = 0.25 — 2.5. 



